Апроксимація сплайнами.
Якщо функція f(х) задана своїми значеннями в (n + 1) – точці, то її можна апроксимувати інтерполяційним поліномом Лагранжа або Ньютона Рn(x) степені n . Але при великих n виникають пульсації многочлена Рn(x) між вузлами (оскільки він має точки максимуму і мінімуму у дійсних нулях похідної). Усунути це можна зменшенням степені полінома, інтерполюючи тільки частину точок (кускова інтерполяція). Тоді для всієї множини точок одержується результуюча “зшита” або “склеєна” крива, яка неперервна на всьому інтервалі, але у вузлах “склеювання” окремих кривих її похідні матимуть розрив. Похибка такого наближення може виявитися значною. Крім того, таке наближення неможливе у тих випадках, коли важливе значення має гладкість апроксимуючої функції.
У ряді випадків інтерполяційний поліном Рn(x) може не наближатись до f(х) при n → ∞ навіть для диференційованої довільне число разів функцій, наприклад, функції
–функція Рунге на відрізку [-1,1]
для якої
Як видно з рисунку, розходження відбувається в інтервалі 0,726…≤│х│< 1.
Побудова полінома невисокої степені для великого числа точок методом найменших квадратів не забезпечує задовільну відповідність кривої заданим точкам (сума квадратів різниць між значеннями многочлена і функції в даних точках хоча і мінімальна, але значна). Підвищення ж степені многочлена приводить до значного ускладнення в обчисленнях.
Останнім часом в теорії апроксимації широкого використання набули сплайни (spline-рейка, пружний брусок), які дозволяють в деякій мірі усунути недоліки класичних методів наближення.
Гладкість і гнучкість сплайнів обумовлює їх широке використання в машинній графіці, САПР, а також в математиці. Сплайни мають ряд переваг:
простота реалізації на ЕОМ;
висока точність апроксимації одночасно і функції і її похідних;
в ряді випадків дозволяють отримати наближення з мінімально можливою похибкою на даному класі функцій (порівняно з іншими методами).
Інтерполяція кубічними сплайнами. Одержали найбільш широке застосування в інженерній практиці. Це обумовлено тим, що поліном 3-го порядку є найпростішою кривою, яка має точку перегину, що забезпечує її хороші інтерполяційні можливості. Кубічні сплайни мають на всьому відрізку апроксимації неперервні похідні до другого порядку. Така гладкість звичайно достатня для більшості практичних задач. Невисока степінь полінома спрощує обчислення і зменшує похибку обчислення.
Візьмемо n + 1 [х0…xn, у0…уn] експериментальних точок, або точок, які задовольняють деяку функцію. При сплайновій інтерполяції (кубічній) виконуються наступні умови:
через кожні дві сусідні точки проводять поліном третьої степені;
для і-го полінома Рі, який проходить через точки з координатами хі-1, уі-1 і хі, уі є справедливою умова інтерполяції
; (1), (2)
Для n = 4 з’єднання інтерполюючих поліномів Р1(х)… Р4(х) буде таким:
Перехід від полінома Рі, що з’єднує точки ( хі-1, уі-1 ), та ( хі, уі ) до полінома Рі+1, котрий проходить через точки ( хі, уі ) та ( хі+1, уі+1 ), повинен бути плавним (без зламів). Щоб виконалась ця умова, значення першої та другої похідної сусідніх поліномів Рі та Рі+1 в їх спільній точці ( хі, уі ) повинні бути рівні
(3)
(4)
Таким чином, для всіх внутрішніх поліномів формалізовані чотири умови, котрі достатні для побудови цих поліномів.
Обидві кінцеві точки ( на кінцях відрізку апроксимації ) мають тільки по одній сусідній. Тому для кінцевих поліномів потрібні додаткові умови. Будемо вважати, що точки, які лежать зовні інтервалу [ x0, xn ], можна апроксимувати прямими; тоді одержимо ще дві додаткові умови
(5)
(6)
Таким чином, маємо шість умов (1) – (6).
Кубічний поліном будемо шукати у вигляді:
(7)
Необхідно виз...